In this paper, we introduce and solve a system of bi-Drygas FUNCTIONAL EQUATIONs \begin{EQUATION}\left\{\begin{aligned} &f(x+y,z)+f(x-y,z)=2f(x,z)+f(y,z)+f(-y,-z)\nonumber\\ &f(x,y+z)+f(x, y-z)=2f(x,y)+f(x,z)+f(-x,-z)\nonumber\end{aligned}\right.\end{EQUATION}for all $x,y,z\in X$. We will also investigate the Hyers-Ulam stability of the system of bi-Drygas FUNCTIONAL EQUATIONs.